Astronomy & Astrophysics manuscript no. oh-great-final2 
February 25, 2013 



©ESO2013 



Letter to the Editor 



High-resolution absorption spectroscopy of thie OH ^113/2 ground 

state line 

H. Wiesemeyer\ R. Gusten\ S. Heyminck\ K. Jacobs^, K.M. Menten\ D.A. Neufeld^, M.A. Requena-Torres\ and J. 

Stutzki^ 



(N 

o 

(N 

;h 

Ph. 
< 
(N 



< 

o 



> 

cn 

o 



X 



' Max-Planck-Institut fiir Radioastronomie, Auf dem Hiigel 69, 53121 Bonn, Germany 

e-mail: hwiese@mpifr.de 
- I. Physikalisches Institut, Universitat zu Koln, Ziilpicher Str. 77, 50937 Koln, Germany 
^ The Johns Hopkins University 3400 North Charles St. Baltimore, MD 21218, USA 



Received ; accepted 



ABSTRACT 



The chemical composition of the interstellar medium is determined by gas phase chemistry, assisted by grain surface reactions, and 
by shock chemistry. The aim of this study is to measure the abundance of the hydroxyl radical (OH) in diffuse spiral arm clouds 
as a contribution to our understanding of the underlying network of chemical reactions. Owing to their high critical density, the 
ground states of light hydrides provide a tool to directly estimate column densities by means of absorption spectroscopy against bright 
background sources. We observed onboard the SOFIA observatory the "113/2, J = 5/2 «— 3/2 2.5 THz line of ground-state OH in the 
diffuse clouds of the Carina-Sagittarius spiral arm. OH column densities in the spiral arm clouds along the sightlines to W49N, W51 
and G34.26-I-0. 15 were found to be of the order of 10'"* cm"^, which corresponds to a fractional abundance of 10"' to 10"**, which is 
comparable to that of H2O. The absorption spectra of both species have similar velocity components, and the ratio of the derived H2O 
to OH column densities ranges from 0.3 to 1.0. In W49N we also detected the corresponding line of '**0H. 
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1. Introduction 



The first molecule discovered in space was the CH radical 
(Dunham 1937). At that time, the first predictions were made 
regarding the abundances of molecules in interstellar space, 
identifying the hydroxyl radical (OH) as a promising candidate 
(Swings & Rosenfeld 1937). It took more than two decades until 
the latter was detected in the ism (Weinreb et al. 1963), in ab- 
sorption against the supernova remnant Cas A, in the F - I - I 
and 2-2 transitions between the hyper-fine structure split A- 
double levels of OH's ^113/2,7 - 3/2 rotational ground state. 
Comprehensive models of the gas phase chemistry of diffuse in- 
terstellar clouds constructed by van Dishoeck & Black (1986) 
revealed the importance of the OH radical in the network of re- 
actions leading to the formation of oxygen-bearing molecules. 
Unfortunately, OH column densities in these objects are difficult 
to determine from radio lines, due to frequently observed devia- 
tions of the underlying level population from LTE (e.g., Neufeld 
et al. 2002). The highest densities occurring in diffuse clouds 
amount to ~10'*cm"'' (Greaves & Williams 1994), while their 
mean density is -10^ cm"^ (Snow & McCall 2006; Cox et al. 
1988 derived upper limits of a few thousand cm"^). 

Storey et al. (1981) first detected the A doublet line from the 
OH ground state (^ris/i/ = 5/2 ^ 3/2) towards Sgr B2 near 
the Galactic centre with the Kuiper Airborne Observatory, but 
their spectral resolution (250 km s"') proved inadequate to sep- 
arate the absorption by the line-of-sight clouds from that occur- 
ring in Sgr B2 itself. Here we report the detection of one dou- 
blet line of this transition (at 2514 GHz) and of its isotopolog 



'^OH (at 2495 GHz) with GREAlQ onboard SOFIA, in absorp- 
tion towards the giant Hii regions W49N and W51 and the ul- 
tracompact Hii region G34. 26-1-0. 15. Both lines are inaccessible 
for Herschel/HIFI. The three observed lines of sight are within 
a 15° wide Galactic longitude interval. Those towards W51 and 
G34. 26-1-0. 15 cross the near side of the Carina-Sagittarius spiral 
arm, with W51 (i = 49!5) at a distance of 5.41 (-1-0.31, -0.28) 
kpc (Sato et al. 2010) and G34.26-h0.15 at ~2 kpc distance 
(cf. measurements of G35. 10-0.74, which has a comparable ra- 
dial velocity, Zhang et al. 2009). The line of sight to W49N 
({ = 43°2) first crosses the near side of the Carina-Sagittarius 
arm, grazes the Crux-Scutum arm, and then again crosses the 
Carina-Sagittarius arm on its far side, where W49N is located at 
adistanceof (11.4± 1 .2) kpc (Gwinn et al. 1992). 

2. Observations, data reduction and analysis 

The observations reported here were performed with the GREAT 
receiver (Heyminck et al. 2012) onboard the SOFIA airborne 
observatory (Young et al. 2012), as part of the basic science 
programme (flights on 2011 July 26 and November 8). On 
the first flight the receiver's M and L2 bands were tuned to 
2514.317 GHz for the -11^2 / = 5/2 ^ 3/2 group of OH hyper- 
fine structure lines (in the lower sideband) and to 1837.817 GHz 
for the ^ni/2 J = 3/2 <— 1/2 lines (in the upper sideband), re- 
spectively. On the second flight, the M band was also tuned to the 
2494.695 GHz frequency of the ^03/2 7 = 5/2^3/2 transition 



' GREAT is a development by the MPI fiir Radioastronomie and 
the KOSMA/Universitat zu Koln, in cooperation with the MPI fiir 
Sonnensystemforschung and the DLR Institut fiir Planetenforschung. 
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Table 1. List of the observed ^113/27 
"*0H lines 



5/2 ^ 3/2 i^OH and Table 2. Results of the OH hne profile fitting. 



Transition 




Frequency [GHz] 


(a) 


Ae [s-']<''' 








OH 


2H3/2, 7 = 5/2^ 


3/2 




F-- 
F - 
F -- 


- 2" 
= 3- 
--2- 


^_ 


2+ 
2+ 

r 


2514.298092 
2514.316386 
2514.353165 




0.0137 
0.1368 
0.1231 








'*'0H, 2H3/2, y = 5/2 «- 


•3/2 




F-- 
F -- 
F -- 


- 2+ 
= 3+ 
= 2+ 


^_ 


2- 
2- 

r 


2494.68092 
2494.69507 
2494.73421 




0.0136 
0.1356 
0.1221 



'"' Varberg & Evenson (1993, CTnns = 30kHz). The frequencies for 
"*0H were derived by isotope scaling. *' Pickett et al. (1998). 



j==A 



Fig. 1. Level diagram (not to scale) for the ^03/2 OH ground and 
first excited state. The observed '^OH transitions are indicated 
by bold arrows and labelled with their corresponding frequency. 
The observed "*0H transitions are drawn as grey aiTOWs. 



of '^OH. Typical DSB receiver temperatures were 4500 K and 
2500 K for the M and L2 bands, respectively. Total power sub- 
traction was performed by chopping with an amplitude of 90" 
at 1 Hz. The raw data were converted from XFFT spectrometer 
(Klein et al. 2012) count rates to forward-beam brightness tem- 
peratures with the module kalibrate (Guan et al. 2012) as part of 
the kosmajioftware observing software, analysing the data from 
the calibration loads and the atmospheric total power and allow- 
ing us to fit both the wet component (typically pwv= 10-20 jum) 
and dry content of the atmospheric emission and thus to deter- 
mine the opacity coiTection (a few xO.l). Further processing of 
the data (conversion to main-beam brightness temperature, with 
a beam efficiency of 0.58, and averaging with l/cr^^ weight- 
ing) was made with the class software. The overall calibration 
uncertainty does not exceed 20%. 

Thanks to the high critical density of the J - 5/2 —> 3/2 
transition (5.1 x lO'^cm""' for a 15 K gas, collision coefficients 
from Dewangan et al. 1987, Einstein coefficients as given in 
Tab. [U, we can safely expect almost all OH to be in its ground 
state at the density of the foreground diffuse clouds along the 
sight-line. This makes the determination of column densities 
much more reliable than those derived from the /1 18 cm line. 
Likewise, the ^ni/2 7=1/2 level (64 K above the ^113/2 J = 5/2 
level) cannot be substantially populated either, because we failed 
to detect absorption in the simultaneously observed ^111/2/ - 



Source 


(a) 
^0 


fwhm"'* 


' max 


^^OH 




[km 


s-'] 




[10'"' cmr^] 








W49N 


2.9 ±1.9 


13.4 ±0.4 


0.9 


8.8 ±4.6 


(12.3 K) 


12.3 ±1.0 


10.1 ±1.1 


10.1 




36.8 ±3.5 


11.8 ±4.4 


1.8 


1.6 ±0.9 




39.1 ± 1.3 


2.8 ±0.5 


2.3 


0.5 ±0.1 




60.5 ± 0.3 


10.4 ±3.9 


2.5 


2.0 ±0.2 


W51 


7.2 ±0.4 


4.6 ± 2.2 


2.9 


0.66 ± 0.08 


(8K) 


59.4 ± 3.3 
68.2 ±1.4 


16.8 ±4.5 
4.6 ± 2.7 


5.2 
1.9 


5.6 ±2.1 


G34.26+0.15 


12.6 ±0.2 


5.7 ±0.8 


2.5 


1.12 ±0.21 


(9K) 


28.3 ±0.5 


6.0 ±1.6 


1.4 


0.67 ± 0.09 




55.4 ±0.4 


12.9 ±1.0 


4.8 


4.73 ± 1.20 



LSR velocity of component. The first and third hyperfine component 
(hfc) are -1-2.2 and -4.4 km s"', respectively, off the second one. 
FWHM of Gaussian absorption profile, deconvolved from hfc split. 
Peak opacity in the strongest hfc. The saturated absorption towards 
W49N (L'isr< 25 km s"') provides only a lower opacity limit. 
Column density per fitted velocity component. 
3/2 <— 1/2 transition. The spectral profile, in absence of emis- 
sion, is thus is given by 



(b) 

(c) 



(d) 



Tmb(v) = r,nb,c exp 



■zz- 

/=1 j=l 



1J,U 



(1) 



where T^b and Tmb.c are the main beam brightness Rayleigh- 
Jeans temperatures of the spectral profile (here as a function of 
velocity) and of the continuum (in single-sideband calibration), 
respectively, and A^vc and A^hfc are the number of velocity compo- 
nents and hyperfine components, respectively. While uncertain- 
ties in the calibration temperature cancel out in the opacity deter- 
mination, any residual offset of variance cr^^^ j^ in the definition 
of the continuum level leads to an additional uncertainty in the 
derived opacity of crm,s,- - cr^^jJTc. The absorption spectra 



suggest o-n 



-0. 1, which is tolerable in view of the substantial 



opacities. A simultaneous least-squares fit to the line profiles of 
all velocity components (Eq.[TJ with the opacity 



ij,i^ 



, / , — ^A^OH wa exp -4 In 2 

7T Anhv-y^ g\,] \ \ 






2\ 



(2) 



yields A^oh, the OH column density per velocity component. 
Here Aui is the fwhm of the Gaussian component /, guj and 
gij are the statistical weights {2F H- 1) of the upper and lower 
level, respectively, of a given hyperfine component j, fo.ij is 
the offset of its velocity from the line-of-sight velocity of the 
source, and wa = 0.5 coiTects for the fact that only one dou- 
blet line was observed (Fig. [TJ. Owing to the relatively large 
number of free parameters (3A^vc), a simulated annealing method 
(Metropolis algorithm) was used in combination with a down- 
hill simplex method (Press et al. 1992). The former assists the 
minimisation process in escaping from a local minimum, and 
the latter improves the efficiency of the convergence. The veloc- 
ity structure in the para-H20 1 1 1 - Ooo spectrum towards W49N 
(Sonnentruckeret al. 2010) suggests five velocity components as 
a strict minimum, while with more components, the procedure 
would start to fit noise features. The results are summarised in 
Tab.|2] The OH column density can be expressed by the relation- 
ship A^oH [cm^^] = 7.8 X lO'^TmaxAyfwhm [km s"'] as a function 
of the opacity in the strongest hyperfine component and of the 
width of the absorption profile of the deconvolved spectrum. The 
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uniqueness of the solution was tested with a Monte Carlo study, 
yielding the standard deviation of each parameter. 

3. Results 

For the spiral arm clouds there is no ambiguity in the fit results 
(Fig. |2]l. Towards the three continuum sources the absorption 
is saturated, i.e., in the (5,20), (40,80) and (35,70) kms"' ve- 
locity intervals for W49N, W51 and G34.26H-015, respectively, 
where the derived main line opacities and column densities are 
to be considered lower limits. For W49N, this caveat is corrob- 
orated by the "*0H absorption profile (Fig.|3). The i**OH/"^OH 
abundance ratio is not expected to be affected by chemical frac- 
tionation (Langer et al. 1984). The synthesised opacity in the 
W49N spectrum peaks at t= 5.7. Assuming for the '^O/'^O ra- 
tio the value in the 4 kpc ring (327 + 32, Wilson & Rood 1994; 
Polehampton et al. 2005 found no evidence of an abundance 
gradient with galactocentiic distance), the estimated opacity in 
our '^OH detection (t~0.2) would require the main line opacity 
to be higher by at least an order of magnitude with respect to 
that estimated by the absorption profile fit. Unfortunately, in the 
30-40 km s"' velocity range the '^OH absorption is affected by 
a telluric ozone feature, and a confirmation by observations of 
sources with a more favourable velocity is planned to definitely 
rule out a baseline lipple. In the unsaturated wings of the absorp- 
tion profile, the sensitivity of the coiTesponding '^OH measure- 
ment is no longer sufficient to estimate the "^OH/OH abundance 
ratio. A two-component fit to the '**OH absorption (Fig.[3]l yields 
a column density of 4 x 10'^ cm"- for the whole absorption fea- 
ture. Although the spiral arm clouds along the sight-line exhibit 
substantial opacities on the order of unity, the absorption is not 
saturated and OH column densities can be derived whose accu- 
racy only depends on the signal-to-noise ratio, the quality of the 
fitted profile, and the assumed continuum level. For a compar- 
ison of our OH column densities with those observed for H2O 
(Sonnentrucker et al. 2010) and inferred for H2 (Godard et al. 
2012), the absorption profiles of the clouds towards W49N and 
W5 1 are integrated within the velocity intervals of Tab. |3] With 
the exception of a velocity interval with an abnormally low wa- 
ter abundance, probably due to a spectral baseline problem, the 
H2O/OH ratios are in the range 0.3-1.0. Plume et al. (2004) have 
determined the H2O/OH ratio by a comparison of submillimeter 
H2O observations with ground-based radio observations of the 
18 cm transitions within the ground rotational state of OH. They 
thereby estimated a H2O/OH ratio of 0.4 at i'isr=68 kms ', in 
agreement with our measurement of 0.6. The ratio of ~0.3 mea- 
sured by Neufeld et al. (2002) towards W51, at wisi— 6 kms"', 
compares to our value of 0.4 in the L'isr=(-l,ll) kms"' in- 
terval. Their OH column density is compatible with our value 
(8 X 10'^ cm"^ and 6.4 x 10'^ cm"^, respectively). Observations 
of the ^ni/2 <— ^113/2 cross-band transitions towards Sgr B2 
(Polehampton et al. 2005) suggest a H2O/OH range of 0.6-1.2. 
Generally, discrepancies between different sets of data may be 
due to NLTE effects in the radio lines, different spectral resolu- 
tions and definitions of velocity components, and uncertainties 
in the definition of the respective continuum levels. 

The chemistry leading to interstellar OH and water has 
been considered in many theoretical studies over the past 
thirty years (e.g. Draine et al. 1983; van Dishoeck & Black 
1986; Hollenbach et al. 2009, and references therein). Three 
main pathways to OH have been identified in diffuse and 
translucent molecular clouds. The first pathway involves an 
ion-molecule chemistry, initiated by the cosmic -ray ioniza- 
tion of H2 or H. The resulting H2^ and H^ ions can lead to 



15 - OH. ^n^/j, J = V2 ^ ^/i- i'D = 2.514 TH 
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Fig. 2. Top: OH absorption towards W49N. The velocities of the 
line-of-sight clouds and of W49N are indicated. The red dashed 
line is a least-squares fit. Middle: Same for the OH absorption 
against W51e4. The relative positions and strengths of the hfc 
splitting are indicated for the v\st-l km s"' component. Bottom: 
OH absorption against G34.26-i-0.15. 




Fig. 3. '^OH absorption (top, grey-shaded), a least-squares, two- 
component fit to it (dashed line) and OH absorption towards 
W49N. The spectra are scaled by the corresponding continuum 
level, to facilitate a comparison. The insert at the top shows a 
telluric ozone feature (as observed in total power), where the 
calibration is more uncertain. 



OH+ through the reaction sequences H2''(H2, H)H3+(0, H2)OH+ 
or H+(0,H)O+(H2,H)OH+. In clouds with a low molecu- 
lar fraction, the resulting OH^ is destroyed primarily by 
dissociation recombination with electrons. In clouds with a 
high molecular fraction, however, OH^ is rapidly converted 



H. Wiesemeyer et al.: High-resolution absorption spectroscopy of the OH ^IIs/t ground state line 

Table 3. OH and H2O abundances and their ratios in line-of-sight clouds towards W49N and W5 1 . Statistical errors are given in 
brackets. 









W49N 


ui,,. intervals [kms '] 




W51 


Wis, intervals [km s ' ] 






(30,37) 


(37,44) 


(44,49) (49,54) 


(54,72) 


(-1,11) 


(ll,16)<'-> 


(43,50) 


NoH 

'''H20 

[OH]/[H2] 

[H20]/[H2] 

[H20]/[0H] 


[10"cm-2] 

[lO'^cm-^] 

[lO^Ocm-^] 

[lO-**] 

[lO-**] 


7.1(0.5) 
2.3(0.1) 
5.9(1.9) 
12.0(4.0) 
3.9(1.3) 
0.32(0.03) 


11.3(6.3) 
6.2(0.6) 
6.6(2.0) 
17(11) 
9.4(3.0) 
0.55(0.31) 


1.2(0.2) 1.5(0.3) 
0.1(0.05) 1.5(0.07) 
0.9(0.03) 3.2(0.1) 
13.2(2.5) 4.5(0.9) 

1.1(0.6) 4.7(0.3) 
0.08(0.04) 1.03(0.20) 


18.1(1.3) 
11.6(0.4) 
22.6(3.7) 
8.0(1.4) 
5.1(0.9) 
0.64(0.05) 


6.4(0.9) 
2.6(0.1) 
4.3(0.1) 
15(2.1) 
6.0(0.3) 
0.40(0.06) 


0.23(0.05) 
0.6(0.03) 

3.8(0.9) 


3.1(0.6) 
2.5(0.1) 
2.6(0.8) 
11.9(4.3) 
9.6(3.0) 
0.81(0.16) 



Sonnentrucker (2010), ""' Godard et al. (2012). ''^' The sensitivity of the OH observations is not sufficient for this velocity interval. 



to 030"^ by a series of two H atom abstraction reactions: 
OH+(H2,H)H20+(H2,H)H30+. The H3O+ ion then undergoes 
dissociative recombination with electrons to form OH or H2O. 
The branching ratio for this process is important in determin- 
ing the resultant OH/H2O ratio and has been studied in two re- 
cent ion storage ring experiments (Jensen et al. 2000; Neau et 
al. 2000): these suggest that ~74% to 83% of dissociative re- 
combinations lead to OH, with almost all the remainder lead- 
ing to H2O (and less than ~1% resulting in the production of 
O). In diffuse or translucent clouds, both neutral molecules are 
destroyed by photodissociation, which - in the case of H2O - 
is an additional formation process for OH. A second and dif- 
ferent pathway may be important in shocks or turbulent dis- 
sipation regions, where elevated gas temperatures can drive a 
series of neutral-neutral reactions with significant energy baiii- 
ers: 0(H2, H)OH(H2, H)H20. Finally, OH and H2O may be pro- 
duced by means of a grain-surface chemistry, in which O nuclei 
are hydrogenated on grain surfaces and subsequently photodes- 
orbed. The relative importance of these three pathways will de- 
termine the exact H2O/OH abundance ratio, but all three predict 
a close relationship between OH and H2O. This relationship is 
supported by the observations reported here, which indicate a 
good correspondence between the OH and H2O absorption fea- 
tures; detailed modelling, which must await a larger sample of 
sight-lines, will be needed to interpret observed variations in the 
H2O/OH ratio. We note, however, that the lower end of the ob- 
served range of H2O/OH ratios (0.3-1.0) is predicted by models 
for turbulent chemistry (Godard et al. 2009). We note also that 
the observed distribution of OH is quite different from that of 
OH^; the latter is believed to arise primarily in material with 
a molecular fraction that is too low to permit the efficient pro- 
duction of H30^, whereas the former will arise in clouds with a 
substantial abundance of H2. 

Future data of the OH ground state transition and the rela- 
tively high precision of the resulting column densities will not 
only allow us to assess the correlation between the abundances 
of OH and H2O, but also to re-calibrate less accurate OH column 
densities derived from decades of radio observations. 
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